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Generation  of  Turbulence  in  Couette  Plow 
between  Excer.tric  Cylinders. 

1 . )  Introductory 


Couette  flow  is  besides  the  Hagen-Poiseuile  flow  the 
fundamental  experiment  for  the  study  of  the  flow  properties 
of  liquids  and  the  study  of  stability  and  transition  to  tur¬ 
bulence  in  ducts.  Although  the  theory  of  stability  was  most 
successful  by  discovering  the  generation  of  cellular  vorti¬ 
ces  [l]  when  the  centrifugal  forces  act  stabilizing  there 
does  not  exist  a  complete  understanding  of  the  transition 
phenomenon  to  turbulence  which  was  observed  [l, 2, 3*^]  when 
the  centrifugal  forces  act  stabilizing.  An  example  for  the 
first  case  is  a  rotating  inner  cylinder  and  the  outer  cylin¬ 
der  at  rest  and  for  the  second  case  the  outer  cylinder  rota¬ 
ting  and  the  inner  cylinder  at  rest.  Only  recently  it  was 
pointed  out  [5]  that  the  observed  transition  could  have  been 
caused  by  vibrations  or  excentricities  which  were  produced 
by  imperfections  of  the  Couette  apparatus  used  in  the  experi¬ 
ments.  The  use  of  ball  bearing,  large  ratios  of  the  length 
to  the  diameter  of  the  cylinders,  large  dimensions  and  cylin¬ 
ders  bent  from  sheet  point  to  the  likelihood  of  such  imperfec¬ 
tions.  When  annihilating  excentricities  and  vibrations  it  in¬ 
deed  was  found  [5]  that  the  flow  is  completely  stable  up  to 
Reynolds  numbers  Jet  not  attained.  Also  the  tncoretical  stabi¬ 
lity  proof  on  the  basis  of  small  two-dimensional  perturbations 
shows  stability  [5*6j*  furthermore  with  determined  excentri¬ 
cities  transition  to  t;  .  .  ulence  w=*c  obtained  with  a  definite 
dependency  on  the  Reynolds  numb  .  There  o_em  to  be  enough 
indications  for  the  ereatio.  of  turbulen .c  by  excentricities. 
It  is  the  purpose  of  this  in .  -stig  .tion  nj  give  a  ‘  ‘-eorctical 
explanation  of  this  transition  phenomenon  ii.  the  ence  of 
excentric  cylinders. 


Transition  can  be  caused  either  oy  flow  instability  or 
by  separation.  In  the  latter  case  there  occurs  near  the  wall 
counter  flow  wich  initiates  transition.  Preliminary  experi¬ 
mental  observations  gave  strong  evidence  for  this  type  of  tran 
sltion.  Therefore  the  Investigation  will  deal  with  the  sepa¬ 
ration  effect  opposite  to  the  original  intention. 

2. )  General  assumptions  and  notations. 

As  mentioned  before  a  rotating  outer  cylinder  and  the 
inner  cylinder  at  rest  will  be  assumed.  Plane  motion  will  be 
regarded.  This  means  that  this  investigation  refers  to  rela¬ 
tively  long  cylinders  so  that  end  effects  will  not  influence 
the  middle  part  of  the  flow. 

The  center  of  the  inner  cylinder  will  be  regarded  as  cen¬ 
ter  of  reference  (fig.  1).  The  excentricity  of  the  outer  cy¬ 
linder  of  radius  r1  will  be  denoted  by  e.  The  ratio  e/ r.  is 
regarded  as  small  so  that  higher  orders  of  this  ratio  can  be 
neglected.  The  discussion  of  the  boundary  conditions  will  show 


that  the  satisfaction  of  these  conditions  at  the  excentric 
boundary  involves  rather  cumbersome  numerical  calculations. 
Therefore  the  outer  boundary  will  be  altered  some  what  in 
the  way  that  the  boundary  conditions  will  be  shifted  from 
the  excentric  to  a  hypothetic  centric  circular  boundary  with 
the  same  radius  as  the  excentric  boundary.  In  the  case  of  a 
small  ratio  of  the  excentricity  to  the  mean  width  of  the  gap 
this  is  of  negligible  influence  whereas  at  larger  excentrici- 
ties  the  boundary  conditions  give  a  periodically  alternating 
in-  and  outflow  at  the  hypothetical  boundary.  One  may  say 
that  these  conditions  are  still  in  the  neighborhood  of  the 
real  boundary  conditions  if  the  excentricity  is  a  larger  frac¬ 
tion  of  tne  average  width  Of  the  gap. 

Navier-Stokes  equations  will  be  considered  with  out  ne¬ 
glecting  any  terms.  However  the  inertia  terms  will  be  lineari¬ 
sed  by  assuming  as  mean  flow  the  Couette  flow  between  centric 
cylinders.  Therefore  only  small  perturbations  should  be  gene¬ 
rated  by  the  excentricity.  This  means  that  the  investigation 
is  restricted  to  excentricitles  which  are  small  in  comparison 
to  the  mean  width  of  the  gap. 

The  linearization  is  not  in  agreement  with  the  physical 
problem  of  separation.  Indeed  separation  Is  an  effect  of  fi¬ 
nite  Inertia  fences  as  fluid  particles  are  subjected  to  a  fi¬ 
nite  deceleration.  Therefore  the  calculations  presented  here 
can  only  show  the  tendency  to  separation.  One  can  not  expect 
that  the  calculated  values  of  the  parameter  which  characteri¬ 
zes  separation  will  be  in  good  agreement  with  experimental 
values.  The  comparison  with  own  experiments  will  show  that 
the  calculated  values  are  to  small.  But  this  lack  in  agree¬ 
ment  characterizes  all  calculations  of  similar  problems  as 
the  mathematical  difficulties  are  invincible  with  out  the  sim¬ 
plification  of  linearization. 

3 . )  Boundary  conditions . 

Tne  circumferential  velocity  of  the  outer  Linder  will 
be  denoted  by  U*.  With  the  center  of  reference  1..  the  center 
of  the  inner  cylinder  one  has  according  to  fig.  2  the  tangen- 


tiai  and  radial  components 

U*t  »  U*  cos  f  U*  .  -  U*  sin  f 

Introducing  the  geometric  relations 

X,  sin  f  •  esin  9  Xicosf*  •  cos  9  -  £ 

one  obtains 

ut  m  u*  $  (  4  -  $  eot<P  ),  uf  -  “  U*  | 

From  (2)  one  derives 

i  •  e  cos  9  ♦  n/  4  - 


*  *» * 


The  series  expansion  is,  if  the  dimensionless  excentricity 


is  introduced 


£  t  t  cos  9  *  i  [  1  -  £  -^1  e*  (l  -cos  2?) J 

With  the  above  mentioned  assumption 

€*  «  4 


(5 

(6 


and  the  notations 


one  obtains 

r  m  f  (4+  A  co 3  ^ 


A  -  £/? 


(7a 


(7b 


This  gives  when  introduced  in  (J>)  the  following  expressions 
for  the  tangential  and  radial  components  of  the  circumferen 
tiul  velocity  U*  of  the  outer  cylinder 

V*  m  U*  ,  U*  »  -  u*  A  sin  <p  (8 


Introducing  the  mean  width 


fc.  *  t  -  n  (9 

and 

*  7*-  *  $  -  4  (lo 

ro 

the  second  equation  (8)  can  also  be  written 


V*  •  -  U*  sin  9  (11 

This  shows  that  excentrici ties  of  the  order  e* €  4 
influence  only  the  radial  velocity  component. 


The  width  of  the  gap  is 

h  «  «j-r,  ♦  €  cos  <p 


(12 


-  6  - 


Then 


h  *  h,  ( 1  «•  ^  cos  9 ) 


03 


Later  tne  variable 


y 


r-  vt 

r. 


(14 


will  be  introduced.  It  is  the  dimensionless  distance  from 
the  inner  wall.  The  value  of y  at  the  outer  boundary  is  8 
according  to  (1o).(13)  gives 

*  S  ^  *  1  *  t  C0S  9  ) 


V.'ith  the  notation 

one  obtains 

8  *  &0(4*<xcos<f) 

As  the  calculations  will  be  restricted  to 

«  •  ft  «  1 

tne  higher  powers  of  8  can  be  approximated  by 

s  £(4 +  2*  cos  <p )  #  6**  8f'\4  +  3MCOS<t)'  ... 


(16 


(17 

(18 


O* 


It  will  be  snown  that  the  limitation  (lb)  is  not  absolu¬ 
tely  necessary.  Nevertheless  it  will  be  Introduced  as  other¬ 
wise  the  calculations  get  to  extensive. 

4 . )  Basic  equations. 

Denoting  by  u  the  tangential  by  v  the  radial  velocity 
(fig.  0  by  p  pressure,  by  y  kinematic  vis-.osity,  by  index  r,  <p 
differentiations  with  respect  to  r,  ?  the  .  vier-f  tokes  equa¬ 
tions  for  the  circumferential  and  radial  direction  refering 
to  the  twodimensional  motion  are 

rvuT  tu«f  *uv  -  ♦  v[uw  ♦  ^  ^  f  A  Vf  ]r  (20 


vvr 


(21 


-  7  - 


Introducing  the  stream  function  y  defined  by 

“  *  -  Vr  .  V  •  -i  VV 

and  differentiating  (2o)  with  respect  to  r,  (21)  with  resj ec 
to  9  then  eliminating  the  pressure  and  introducing  Laplacf s 
operator 

Ay  *  V-r  *  7  H  VW 

one  finally  obtains 

-  V„  A  V*  •  vrAAy  ^22 

Tlie  stream  function  will  be  composed  of  two  parts  y*  of  the 
mean  flow  and  y'  of  the  perturbation 

y  *  y*  *  y' 

As  mentioned  before  for  the  mean  flow  the  Couoote  flow 
between  centric  cylinders  will  u-  introduced  the  velo¬ 
cities  of  which  are 


With  this  (22)  is  transformed  to 


U  A  w'  •  V*'  -  WflW  *  V  r  dd  y ' 


V  »0 


(24 


Linearizing  tnis  equation  by  neglecting  the  second  order  ter 
in  y‘  writing  y  for  y’  and  introducing  trio  abbreviation 


one  obtains 

K(£i-1  )dvV  *  r4Ady  (26 


The  variable  y  defined  by  (14)  wiii  be  In.  reduced.  Tnen  .  .e 
equation 

K(zyy*)[«* y)VyyV  ♦  M*yHw*ym J «  [(^y)*v^ry  ♦  “ 

U, 

-w*y»*Wr  *"*y>v,*+Vf9  *w+w«y»Stm  ♦Vw,J 


1  o  aKV  •»  4  nrt  H 


b . )  Solution  for  zero  Inertia  terms. 


According  to  the  geometry  of  the  problem  a  solution 
,riodic  in  the  circumferential  angle  y (f ip. * )  is  to  be 
expecteu.  Solutions  with  this  property  are  well  known  for 
vanishing  inertia  terms  as  the  basic  equations  then  are 
reduced  to  the  biharmonic  differential  equations.  One  has 
with  tno  dimensionless  radius  r  «*  r/rQ 

y  »  a4  ♦  r  *a,rl  *ccs  y  (b4r  ♦  b \r4+  btrl  *  b,r  Inr) 

+tasl<f(c4vx*  c,*1 1  ctr*  *■  c*  ) 

+co$3y (d<r,id4r*1  ♦d,r'  ♦  d^r’4  ) 


.•’roni  this  one  obtains  the  velocity  components 


u  -  -  yr 

*  -  |  axr'*  ♦  Z  a,r  *  cos  qp  ( b*  - 4^ r*1 .  3b* r*  * 6^  ♦  fe, /nr  ) 

♦  cos  2y  (c4-Zr  -  2cx v*  ♦  bcjr*  ) 

♦  co s3f  (3d,r*  -3d,Y'*  * 5  d}r*  -  d,r** ) 

....  j 

v  -  y  V v  *  -  {  *i*  y  (fc«  ♦  ♦  b,r*  +  b,  In  r  ) 

♦  CStn  i<p(c4r  *C 4r*f  W  C,r:  ♦  %?"* ) 


(28 


♦  3$<n  If  (<t4r*  *d/'f  *d,  rf  +  d+r'*) 

..  , 


(29 


The  boundary  conditions  u»  v  «  0  at  v«  4  give 

*>4  ♦  K  ♦  ^  •  0  a,  ♦  2  aj  •  0 

c,  *  c*  .  c,  *  c,  -  0  b4  -  b*  ♦  3  ^  ♦  b,  •  0 

2c,  -  2c*  ♦  *€,  •  0 


Introducing  (7#B)  in  (28,29)  one  obtains  for  the  outer  boun¬ 
dary,  when  powers  of  er  are  developed  into  power  series  of  A 
and  only  first  order  terms  in  A  are  regarded 

-  U*  •  a*  q4(*  -  A  co*  9)  ♦  2  a,  9  (w  *  A  cos  «f  ) 

♦  cos  v  [b4  *  b*^1^  -2  Acos  p)  ♦  2  Acoj  ) 

♦  b^faWnf  )  |n(4  fAcosyJ 

♦  co$2f  [20,9(44  Acos? )  -  2c*?1 0-3  A  cos  <p) 

♦  *c,9*(4  4  J  Acos  ?)J 

♦  *  « 

U*Asi«?  ■  si*9[b<4.|\9i(‘f-lAco5f )  ♦  b,?  (1*2Xcos<p)  *b*  Jnf 

♦  b,  ln(  -f  ♦  Acoj  ?)  J 

♦2sm2<p[c,f  (4*Ace s  ?)  «■  c*  9*J  (4-3A  cos  ?) 

4  C9(f(4+3  A  cos  v )  ♦  c,  $'•*(*  -  A  cosy)] 

Rearranging  terms  with  respect  to  multiples  of  ?  one  obtains 

-U*  -  0*9*  ♦  2a,?  ♦  b,^r*A  *3b,f'A  *f  b,A 

♦  cos  ?  [-a4f**A  ♦  2a, ^ A  ♦  fc,  -  b,?*4  *  3b, 9*  4  b,(1*  I*  ?) 

♦  C.9  A  ♦  3c*  f'*A  ♦  6  c,  q*\] 

♦  cos  *  3fc,  9*A*  £  Ab¥  4  20,9  -  2c*f*J  4  *0,9*4  •••  J 

♦  •  t 

V*Aj*«*  -  s»'*»?[b,*  b*?*4  4  4^1*9  4C,9A  -3c*9*»A  ♦3clAf,-c,A^M] 

♦  stHiyf-b*^  ♦  b,9*A  *  |b,A  *2c,9  4  2c, 9*  ♦  2c, 9*4  20,9%  ...  J 


Equalizing  terms  with  the  same  multiples  in  9  on  the 
right  and  left  side  one  obtains  5  equations  which  together 
with  the  5  equations  (30)  for  the  inner  boundary  determine 
the  1o  constants  ag,  a^»h^,b^>  b^,  c ^ ,  Cg,  c^. 

The  numerical  evaluation  of  the  constants  is  given  in 
table  I  for  A  **  0,08.  Terms  including  6y  are  considered  to 
show  the  convergence.  Cne  sees  that  it  is  sufficient  to 
consider  coefficients  up  to  d  that  means  up  to  the  terms 
containing  3<p  . 


This  solution  surely  would  not  be  sufficient  to  show  the 
separation  effect.  It  merely  should  demonstrate  the  satisfac¬ 
tion  of  the  conditions  at  an  excentric  boundary.  However  this 
solution  is  part  of  the  solution  whicn  considers  the  lineari¬ 
zed  inertia  terms.  This  will  be  shown  later. 


a2~  ',792o1125p 
a.,=  -0,896005614 
bj«  -2,1o384l4oo 
h2~  7# 431b28o>ij 
V  -5027986714 

b^=  25,519629761 

c 1=  -5,0^7627428 
c2=  -1,135904752 
c,=  0,955861533 
w  «  3,22767o0:-2 


Table  I 

d^  o,38l917o28 
dg=  o, 167656633 
<ly  -0,19872580 9 
V  -0,350847855 

e,-  -0,023912624 
e^=  o,o2242o572 
ay  0,023559611 
e,,=  -o,o22o47551 


f  -0,019382031 
fg=  -0,036423168 
fy  o,oo822o991 
f4»  o,o475842o8 
gy  0,011262653 
Gg=  0,018200651 
Zy  -0,006352103 
k4=»  -0,023111201 


6 . )  Solution  considering  linear i z ed  in e r  t i a  terras. 


Similar  to 
expression 


solution  for  zero  inertia  terms  the 


* 

V  *  £  i  (y)  cos  *9 


t 


9$  (y)  s*n  s? 


will  be  introduced  with  the  power  series 


i  (y)  *  £  a*  y*  ,  gt  <y)  •  £  m*  y*  (32 

The  coordinate  y  i3  defined  by  (14).  Introducing  this  expres¬ 
sion  into  the  basic  equation  (27),  then  putting  the  terms  of 
cos  gf  ,  sin  up  to  zero  and  comparing  equal  powers  of  y  one 
obtains  equations  with  which  the  coefficients  of  each  series 
can  be  expressed  by  the  first  four  ones.  To  satisfy  the  con¬ 
ditions  at  the  outer  boundary  one  has  to  introduce  the  expres¬ 
sion  (17#  19)  for  y  and  its  powers.  Powers  and  products  of 
cos  tft  sin  $f  will  occur  which  can  be  expressed  by  the  sinus 
and  cosine;  of  multiples  of  the  angle  <p  as  shown  in  section  5 
for  zero  Inertia  terms.  Then  each  s-  term  in  (32)  would  demand 
the  determination  of  four  constants  from  the  boundary  condi¬ 
tions  . 

These  evaluations  Indeed  would  be  exceedingly  cumbersome. 
Therefore  the  limitation  (lb)  will  be  introduced  so  that  the 
excentrlc  boundary  may  be  replaced  by  a  centric  boundary  with 
radius  r^.  Now  instead  of  (31)  the  expression 

V  «  e  [f  <y>  COS  <p  ♦  g  (y)  gin  <p  ] 
with 

f  •  Z  «.  y‘  .  9  •  I>.  y‘ 

•  • 

is  sufficient. 

The  following  procedure  is  similar  to  tne  one  mentioned  be¬ 
fore:  Putting  the  terms  of  cosy  and  of  siny  to  zero  and 
comparing  equal  powers  of  y  one  obtains  equations  which 
allow  to  express  the  aeries  coefficients  by  the  first  four 
ones  of  each  series.  Tnis  is  a  similar  method  used  by  Oort- 
ler  to  treat  the  free  boundary  layer  flow  neai  a  corruga¬ 
ted  plate  (73  . 

The  following  expressions  are  obtained  for  the  various 
coefficients  a,  m.  The  laborious  derivation  will  be 

here. 


(33 

(34 


omitted 


Following  this  way  one  first  obtains  by  introducing 
(33)  in  (27)  the  differential  equation 

{co*?  -  3f‘  *3f) 

♦  iin  ?  i-g"  -  *3^*  -  3 g'  *  3g)J 

♦  y  lcos<f[-¥f~-6fm+  6f-3f'  +  K(-Zg  .2g’-2g')] 

♦  Jin?  C-*^-  69-  ♦Sy'-V  *K(lf  -Zf  +  in]} 

♦  y*  {  cos  f  f-  6^-67*  ♦  3/"  ♦  K  (-9  ♦  39*  ♦*** >/ 

♦«n?r-69-'-^  ♦  K  (  f-Sf-Sf)]} 

*y’(cosvC-*r-ir*K(g,+  *9')] 

+*r[-¥g*  -29~  .K(-f'-*f'))} 

♦y*  {coj  ?  f  -  {"  *  Kg']  +  tin  <f[-q"  -Kf*  1  )  •  0 


Then  introducing  the  series  expansions  (34 )  for  f,g  and 
expressing  the  series  coefficients  by  the  first  four  ones 
of  each  series  the;  following  expressions  are  obtained,  when 
the  conditions  u  =  v  -  0  .xt  the  inner  boundary  are  introduced 
(See  (44)). 


+  o»25 


(35 


a 


5 


in 


o,45 


a  7 


ap  + 

0,23  d  +  Ki_  0,0;  -  ) 

m2  2 


(36 


a6 

ra6 


ci  -j  p  1  ^  7 

-o,423  r  +  o,2p  +  K (-  0,03 
'3 


+  0,016 


in  ^ 

4J 


(37 


av 


-f-f 


o,4lo7l4  J  -  0,25  2  +  K(+  0,0380952  J  -  o,ol42857m2) 

3  2  *3  a2  (j| 


®  -  -0,4017857  J>_  +  0,25  2  +  K(-  0,0392837  ^  +  o,u1>>52^2) 

*■  3  2 

-  K2- 0,00079363  a2  - 

;n2 


"  °*39583  5  -  0,25  2  ♦  K(  +0  ,0392857  3  t  o,o;24oo79m2 

-'3  2  3  u2 


2  ^ 

+  K  (-0,000661375  :i3  +  o,oo1u5o2  2  ' 

rrl  ' 


-  -0,3916  3  +  0,23  2  +  K  (-  0,0; 80690  3  +  0,0119544  m'- 

' j  2  "  x . . 

+  K  (0,00128307  0,00 1190476  2)  (4. 
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These  expressions  show  clearly  the  influence  of  inertia 
forces  which  is  represented  by  the  terms  dependent  on  K. 
Indeed  in  the  expressions  for  the  coefficients  the  first  two 
terras  which  are  independent  of  K  represent  exactly  the  be¬ 
fore  mentioned  solution  for  zero  inertia  terms.  One  confirms 
this  easily  by  replacing  the  variable  r  by  y  and  by  develo¬ 
ping  in  power  series  of  y  .  One  notices  that  the  series  are 
not  absolutely  convergent  in  K  and  this  means  according  to 
the  definition  of  K  also  in  the  Reynoldsnumber.  The  first 
influence  of  inertia  forces  occurs  in  a^,  ag,  ra^,  m^  with  the 
first  power  of  K  .  Then  in  the  three  following  expressions 
for  the  pairs  a,  m  the  quadratic  power  of  K  is  added.  The 
next  three  expressions  contain  also  the  third  power  and  so 
on.  Therefore  this  series  expansion  only  can  be  used  for  re¬ 
latively  small  R  ynolds  numbers.  This  is  a  well  known  pecu- 
larity  which  for  example  also  occurs  when  applying  the  method 
of  successive  approximation  f8]  . 

The  expressions  (35  to  43)  refer  to  the  boundary  condition 
u  =  v  =  0  at  the  inner  cylinder  where  y  *  0  .  This  gives 

f(o )  •  f'(o)  •  g(o )  -  g'(o)  »  0  (44 

what  means,  that 

a,  m  a«  ■  m#  •  >**,  •  0  (45 

It  may  be  added,  that  for  arbitrary  boundary  conditions 
at  y  =  0  the  expressions  (35  to  43)  contain  the  following 
additional  terms 
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The  first  boundary  condition  in  (S)  for  the  outer  boundary 
y  mi  is  satisfied  by  the  mean  flow  whilst  the  second  one  in 
(8)  has  to  be  satisfied  by  the  secondary  flow.  Hence  the  boun 
dary  conditions  for  the  secondary  flow  are 

,  .  .  (46 

u  ■  0  t  v  •  -  U*  ■■ 1 6  ■  -  sin  <p 

6  +  1 

Introducing  (33,34)  one  obtains 

(47 

f'Ot )  •9,( <f)  »0  ,  9 f<T)  -  0  ,  .  U* 

The  boundary  conditions  (44,47)  give  eight  equations  to 
determine  the  eight  constants  aQ,  a^,  a2,  a^,  mQ,  m^,  m2,  rn., . 
Ar  mentioned  above  the  conditions  (44)  yield  a  *  a.j  »  inQ=*  :rij 
So  that  the  four  equations  (47)  still  have  to  be  solved.  This 
work  was  carried  through  numerically  with  the  aid  of  computer 


7. )  Numerical  calculations. 

The  numerical  quantities  Inserted  for  4*9  (1o)  are 

,  f  •  4.1 

Due  to  the  before  mentioned  semiconvergence  of  tne  seri 
expansion  there  exist  certain  limits  for  the  Reynolds  numb 
beyond  of  which  the  convergence  is  not  any  more  satisfactory. 
With  the  assumed  value  4=  o,2  the  limit  is  K  *=  1o^.  With  (2t) 
this  corresponds  to  a  Reynolds  number 

ftt  •  »  *  40* 

The  calculations  were  performed  f'  r  two  values  of  K,  espec- 
t.ively  Re.  The  boundary  conditions  (47)  give  the  it  llo.’  z 
expressions 

1.  )  K  *  3  •  1o>.  Re  =  1,32  • 


f(8)  m  0,016011  ai  *  0,005206  a,  ♦  0,025806  +  0,00*302  m,  *  U* 

$(<f)  «  -0,025/06  a,  -  0,00*302  a,  «■  0,026011  mt  +  0,005206  m,  *  0 

f'(8)  --0,15057  a,  ♦  0t015S7  a,  ♦  0,5*3 *5  w»,  ♦  0,097*5  m,  «  0 

g'(<0»  -0, 5*3*5  a*  -  0,097*5  a,  -  0,15057  m,  ♦  0,01557  m,  •  0 

The  solution  is 

a_  *  56,09833  U*  su  =  -54,48119  V* 

2  2  (48 

«  -161,82 191  U*  r.3  =  416,2653°  U* 

2. )  K  -  10**  Re  «  4,4  •  105 

« 

f (<f)  -  -  0-MII05  a,  -0016152  a,  -  0, 030515  m,  -  0  0155*1  m,  •  U* 

9(«f)«  q 030515  a,  ♦  0,0155*1  a,  -  0,111905  m,  -  0,016192  m,  «  0 

0,579207  a,  *  0,090829  at  +  0,*13*92  mt  +  0,1*6901  wi3  -  0 

g'(<J)«-  0<H3m  a,  -  0,1*6902  a,  +  0,579207  w,  ♦  0,090929  m3  -  0 

The  solution  is 

a2  *=  -3o,4o^o9  U*  m2  =  2,63489  U*  (4^ 

.  119,69097  U*  =  38,27269  U* 

Now  the  excentricity  can  be  evaluated  for  /hieh  sepa. 
tion  will  occur.  It  is  to  be  expected  that  separation  first 
occurs  at  the  inner  boundary  y  =  0  .  The  beginning  of  separa¬ 
tion  will  be  characterized  by  the  zero  value  of  the 
derivative  of  the  total  circumferential  velocity  In  radial 
direction.  This  means 


U'(o)  ♦  u  (o)  -  0 


(5° 


(23)  gives 


u<>’  -jfcrr  u*  <»•«-  ?£-> 

*j*Mr  u* 

Inserting  &  -  o, 2  one  obtains 


U'(y)  -  2.7Z„  U*  (-»  •  pSy|>) 
l/*(o)  -  5,4?...  U* 

By  differentiation  of  (33)  one  derives 


(31a 


(31b 


-  -£{  f*<y>  cosy  ♦  $*(y)  sin  y} 

«  -  t{f2a,  ♦>  6a, y  ♦  ...  J  cos  y  ♦  [2m,  *  6m,  y  *  -  ]**y/ 
u(0)  •  -  t {  2 a*  cos  y  ♦  2  m,  sin  y  } 

(52 

Introducing  (51,52)  in  (5o )  one  obtains 

5,45  U*  -  £  (2 a,  cos  y  ♦  2m,  sin  y  )  *  0 

The  coefficients  were  calculated  before  for  two  values  of 
K  in  (47,48).  Therefore  (53)  determines  the  critical  excen- 
tricity  for  certain  angles  y .  The  smallest  excentriclty  ob¬ 
tained  should  be  regarded  as  the  critic  .1  value.  It  is  suffi¬ 
cient  to  calculate  the  excentricity  in  eec  .  case  for  two  an¬ 
gles  . 

1 .  )  K  «  3  •  lo-5 

y  -  0°  5,45  Uf  -  2  £  a2  =  0 

t  =  0,04867,  e  -  1,o  13^2  in.n 
y  ••=  270°  5,45  U*  +  2  im,  =0 

£  -  0,05006^  e  »  1,o4Z"  o  ;nm 

2. )  K  *  104 

y  =  180°  £  =  0,  o89ri  w  ,  e  -  1,86873  mm 

y  =  90°  £  -  1,03506  .  e  -21,36375  mm 


8. )  Comparison  with  experiments . 


To  prove  the  theory  experiments  were  performed 
with  a  Couette  aparatus  the  cylinders  of  which  were 
excentric.  The  dimexisioxis  of  tne  apparatus  were  the 
following 

rQ  =  21  aim  ^  =  1 , 1  y 

r^  =  23  *mn  6  *»  o,  19 

e  =  o,t>;  1;  2;  2,t>i  3  .nm  e  =  o,o4y6  (e  =  1mm 

As  shown  in  fig  3  no  ball  bearings  were  used  for  the 

support  of  the  inner  cylinder.  The  occui’ance  of  separation 
was  observed  visualy  w.th  a  technique  described  earlier  1^1  . 
With  the  excentricity  e  =  1,p  mm  there  was  found  a  critical 
Reyxiolds  number  for  separation  Re  =  U*»*/v  =  1,0:*  .  1o^ 

[  ‘j  I  ,  These  experiments  were  now  extended  to  the  before 
mentioned  excentricities .  With  e  =  o,p  mm  no  separation 

C 

could  be  observed  up  to  the  Reynolds  number  4  .  lo  .  The 
critical  Reynolds  numbers  whxch  were  observed  are  plotted 
in  fig  4.  As  thxS  figure  shows  there  exists  a  definite  de¬ 
pendence  on  the  excentr.'.city .  This  confirms  the  supposition 
that  turbulence  is  generated  by  separatioxi. 

The  calculation  fox-  Re  =  1,3  •  1o^  gave  separation  with 
fc  =  u,o48.  This  exceutx'icity  was  realized  m  the  experi. neats 
and  as  fig  4  shows  the  corresponding  Reyxiolds  number  for 

4 

separation  is  2, 1  .  lo  .  One  sees  that  the  theoretical  value 
is  16  times  to  small.  This  seems  to  be  a  satisfactory  abi>ee- 
•lexit  for  a  first  oi'der  approximation. 

It  may  be  mentioned  that  earlier  calculations  of  the  se¬ 
paration  at  a  corrugated  plate  showed  a  sensitive  influence 
of  the  corrugation  1 7  ]  •  Now  accordixig  to  the  experimental 
comparison  this  sexisitivity  seems  to  be  more  ixu'luenced  by 
the  degree  of  approximation  than  by  physical  effects.  This 
..lakes  an  evaluation  of  tne  inertia  forces  necessary.  * 
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